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We have applied the Green-Kubo formula to investigate the shear viscosity in the SU(2) light- 
flavor quark matter at finite temperature under the external strong magnetic field e\B\ ~ m 2 . 
For this purpose, we employ the temperature-modified instanton model and the Schwinger method 
to induce the magnetic field. The quark spectral function with the finite width motivated by the 
instanton model is adopted to compute the shear viscosity. We find that shear viscosity increases 
as temperature increases even beyond the transition temperature To = 170 MeV if temperature- 
dependent (TDP) model parameters is used. On the other hand, with temperature-independent 
ones the shear viscosity starts to drop when temperature goes beyond To. Although the magnetic 
field reduces the shear viscosity in terms of the magnetic catalysis, its effects are almost negligible in 
the chiral-restored phase even for very strong magnetic field, e\B\ ~ 10 20 gauss. We also compute 
the ratio of the shear viscosity and entropy density n/s. Our numerical results are well compatible 
with other theoretical results for a wide temperature regions. We obtain the parameterization of the 
temperature-dependent ratio from our numerical result as n/s = 0.27 — 0.87/t + 1.19/t 2 — 0.28/t 3 
with t = T/T for T = (100 ~ 350) MeV and e\B\ = 0. 
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I. INTRODUCTION 



With the energetic developments of the heavy-ion collision (HIC) experiments by Relativistic Heavy-Ion Collision 
(RHIC) at BNL and Large Hadron Collider (LHC) at CERN, the properties of quark-gluon plasma (QGP) have been 
widely investigated. One of the most highlighted observations from the HIC experiments is that QGP behaves as 

i-^ almost perfect fluid, being characterized by the small lower-bound value for its shear viscosity, i.e. the Kovtun-Son- 
Starinets (KSS) bound [4]: f > j^- It has also been supported by theories, such as the viscous hydrodynamics p] 
and AdS/QCD [HI3j. This observation implies that QGP is a strongly coupled system [5j. The viscous hydrodynamic 
simulation for the elliptic flow V2 with the Monte-Carlo (MC)-Glauber initial condition reproduced the Au+Au collision 

■^- data with 2 — i j while the value 2 w ^ did the experimental data with the MC-KLN initial condition [TJ [5]. 

These observations indicate that the initial condition for the hydrodynamic simulations is a potential uncertainty 
to determine the shear viscosity. It is also worth noting that all the hydrodynamic simulations have employed the 
temperature-independent shear-viscosity values so far. Hence, the information of the temperature dependency of the 
^. shear viscosity is important to estimate the QGP shear viscosity. In addition, the initial quantum fluctuations, such 
as the color-charge fluctuation, can also cause the uncertainty as well [7]. For a recent status for the shear viscosity, 
one may refer to Refs. 8,9. 



The shear viscosity can be theoretically explored by the Green-Kubo formula in terms of the linear response 



theory 10 20 . Since QGP is a strongly coupled system, it can only be studied via nonperturbative methods in 
principle, such as low-energy effective QCD-like models or lattice QCD (LQCD) simulations. From the effective 
models, such as the Nambu-Jona-Lasinio (NJL) model, the shear viscosity has been scrutinized extensively as a 
function of temperature (T) and/or quark chemical potential (/z) [TTJ [TTJ [T!5] . In Refs. [HIE?], it was argued that the 
quark spectral function in a simple local mean-field (MF) only gives trivial results, i.e. 7/=0. To overcome this difficulty, 
one needs to go beyond MF. One of the remedies for this issue is to consider a finite width for the quark spectral 
function [TT]. The shear viscosity was also investigated by LQCD simulations 21j, dissipative hydrodynamics [TOlllSj . 
chiral perturbation theory (xPT) [T2], perturbative QCD (pQCD) [13], and holographic models [2"1 131 12"2"] . 

In addition to the shear viscosity of QGP, the effects from the external magnetic field produced in the peripheral 
HIC experiments have also attracted much attention [33J . Although the produced magnetic field is reduced by a factor 
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~ 10 4 after ~ 3 fm/c [53], its strength is still very strong in the order of pion mass squared: e\B\ ex m\ ~ 10 18 gauss. 
Such a strength is comparable to the magnetic field of neutron stars. Moreover, recently people have speculated that 
a strong external magnetic field may generate the chiral magnetic effect and chiral magnetic wave which will induce 
in the CP-odd domains in QGP [35] • Moreover, the magnetic field to the QCD matter enhances the spontaneous 
break-down of chiral symmetry (SB^S), i.e. the magnetic catalysis |26) . Hence it is interesting to understand how 
the external magnetic field will change the shear viscosity of QGP. 

In the present work, we want to investigate the shear viscosity of the SU(2) light-flavor quark matter at finite 
temperature under the strong magnetic field. For this purpose, we employ the dilute instanton liquid model (LIM) for 
the light flavor SU(2) sector [27J I2E]- This model manifests the nontrivial quark-instanton interactions via the quark 
zero mode, resulting in the natural UV regulator by construction. Since we are interested in the system at finite tem- 
perature, we modified the LIM parameters, such as the average inter- (anti)instanton distance (R) and (anti)instanton 
size (p), using the caloron solution for the Yang-Mills equation with trivial holonomy, i.e. the Harrington-Shepard 
caloron |291 I30j . Although we do not have deconfinement order parameters in our present theoretical framework, 
we consider that the chirally-restored phase can reflect the information of QGP, since the chiral restoration and 
deconfinement temperatures are almost the same for the system with the finite current-quark mass [31H33J . Since 
the quark chemical potential is expected to be small inside QGP created in the HIC experiments, we choose p=0 
throughout the present work for brevity. As mentioned above, the Green-Kubo formula is employed to compute the 
shear viscosity in terms of a quark spectral function |llj . We construct a quark spectral function with a finite width 
A ~ 1/p motivated by the instanton physics. The external magnetic field is induced to the system in hand by the 
Schwinger method |3"3H3"B] . 

Before studying the shear viscosity itself as a function of T, we have to investigate the chiral phase transition within 
the present theoretical framework since the constituent quark mass as a chiral order parameter plays a crucial role in 
the temperature dependency of the shear viscosity. By computing R and p as functions of T with the caloron solution, 
we observe that they are all decrease with respect to temperature, signaling the (partial) restoration of SB^S [3T)ll3"T] . 
At T m Tq, where Tq indicate the chiral phase transition temperature, there appear about 10% decreases in R and p 
in comparison to their values at zero temperature. Using these results and the thermodynamic potential of LIM in the 
leading 1/N C contributions, we explore the chiral phase transition and shows the second order (To sw f 66 McV) and 
crossover (To ~ 170 MeV) transitions for the chiral limit and finite current-quark mass, respectively. It is consistent 
with the universal restoration pattern. 

The numerical results for the shear viscosity are given as functions of temperature as well as the strength of the 
external magnetic field with the temperature- dependent parameters, p{T) and R(T) (TDP), and the temperature- 
independent parameters, p(0) and -R(O) (TIP). With TIP, the curves for the shear viscosity increase up to To then 
decrease smoothly. On the contrary, the curves keep increasing for TDP beyond To, signaling the reduction of 
nonperturbative effects. We also observe a tendency that the external magnetic field lowers the shear viscosity, due to 
the enhancement of the SB^S, i.e. the magnetic catalysis via the magnetic field. According to the second-order chiral 
phase transition in the chiral limit, there is no differences in the shear viscosity due to the magnetic field beyond To, 
since the magnetic-field effects are proportional to the constituent-quark mass in the present theoretical framework 
which is a chiral order parameter. If we go beyond the chiral limit, manifesting the crossover chiral phase transition, 
we observe finite differences in the shear viscosity curves even for T > Tq, but they vanish gradually from T « 220 
MeV. In general, the effects from the magnetic field to the shear viscosity are less than 10% for e\B\ < 100 x to 2 . « 10 20 
gauss. 

We also present the ratio of the shear viscosity and the entropy density 77/s as a function of temperature and the 
strength of the external magnetic field. We find that 77/s decreases smoothly and approach the KSS bound for TDP, 
whereas TIP result undershoots the bound. Moreover, the effects from the magnetic field becomes almost negligible 
beyond T , although the effects are still visible below T . We also compare the present numerical results for n/s with 
other theoretical estimations from the NJL model, LQCD, and xPT, resulting in qualitatively good agreement with 
them. Typical values for the shear viscosity at T = To are given as 1} = 0.02 GeV 3 and 77/ ' s = 0.29 from the present 
model for T = (100 ~ 350) MeV and n B = 0. 

The present work is organized as follows: In Section II, we provide theoretical framework for computing the 
shear viscosity of quark matter: the Green-Kubo formula, finite-width quark spectral function, Schwinger method, 
thermodynamic potential, and so on. The numerical results and relevant discussions are given in Section III and the 
final Section is devoted for the summary and conclusion. 

II. THEORETICAL FORMALISM 

In this section we briefly introduce our theoretical framework. 



A. Shear viscosity at finite temperature 



First, the static shear viscosity 77 is defined according to the Green-Kubo formula, [11] 
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(1) 



where 11^ stands for the retarded (R) quark correlation function and u> for the frequency of the system. The retarded 
correlation function is related to the correlation function as follows: 



nl(u) = YLl(iw 



\iw— >uj-\-ie • 



(2) 



Here w is the fermionic Matsubara (M) frequency. 11^ is defined as the time-ordered tensor current correlator, 



n M(* w ) = - 



1/T 



dr{0\T[J xy (r,T),J X y(0,0)} |0), J xy = - [$fr y d x il>) - (d x ^) lv ^] , (3) 



where r and ip stand for the Euclidean time and quark field. T denotes temperature. One can evaluate 11^ with the 
full quark propagator S by using the fermionic Matsubara formula with w n = (2ri + 1)ttT: 

f d 3 k °° 
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2iri J (2tt) 3 



n F (z)Tr c j^ [k x -f 2 S{iw + z, k)k x ~/ 2 S(z, k)} . 



(4) 



The trace runs over the color (c), flavor, (/) and Lorentz (7) indices. From the first to the second lines in Eq. Q, we 
have employed the fact that the poles of the Fermi-Dirac distribution: 



n F (z) = 



1 



are located at z~i(2n + 1)ttT. The spectral function is related to the quark propagator in the following way, 

s(k ,k ) =rp p p k i. 

J-oo 2?r fc - w 
Hence we obtain the following expression of the shear viscosity in terms of the quark spectral function p{k): 
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Here n' F is defined as, 
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If we adopt the spectral function associated with the free quark propagator with a constant quark mass m given in 
Ref. [37]: 



p(w, k) = 27rsgn[iy] (70^ — 7 • k + m)S(w k — m ). 



(9) 



(This quark spectral function satisfies the normalization condition ^- J p(w,k)dw = 70 [llj as shown in Appendix.) 
Due to the 5-function in the spectral function in Eq. Q, the shear viscosity for the free quark at the mean-field level 
becomes zero, i.e. lim e _>.o J dwf(w)S(w + e)S(w) — as long as f(w) is a regular function. To overcome this problem, 
we introduce a finite width for the quark spectral function as Refs. PHUT]. Thus, we replace the delta function in 
Eq. J9| with a Gaussian functions with a finite width: 



S(w 2 -k 2 - to 2 ) = S{w 2 - E 2 ) 
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T(w,k), (10) 



where the energy for a quark and the momentum-dependent effective quark mass are defined as in the previous 
works [38] by 



Eu 



k 2 + M 2 , 



M k = M (T) 



2 + p 2 (T)k 2 



(11) 



Note that the A ~ 1/p was taken as the width for the Gaussian function. Since the p is only an interaction range 
parameter in this model, it is natural to identify it as the finite width for the quark spectral function. 

It is worth mentioning that M k presents the nonlocal (momentum-dependent) interaction of the quarks. The 
parameter n in Eq. (11) will be determined in such a way that pfw reproduce a typical physical quantity such as 



the chiral condensate. Note that the constituent-quark mass at zero virtuality Mq and average (anti)instanton size p 
are functions of temperature here and will be discussed in detail below. Merging Eqs. (10) and Q, we arrive at the 
following expression for the finite- width (FW) quark spectral function: 



PFw(w,k) = 27rsgn[w;] (j w - 7 ■ k + M k )F(w,k). 



(12) 



Note that the current quark mass m has been replaced by the momentum-dependent effective quark mass asm-) M k 
to regulate the quark-loop integral, in which M k denotes M k + m. pfw a ls° satisfies the normalization condition for 
the quark spectral function as shown in Appendix. The chiral condensate can be related to the spectral function in 
Minkowski space: 
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Performing the Wick rotation for the temporal direction and integrating over (w, iko, k), one is led to 



(2^)4 



1,2 



(13) 



(14) 



To reproduce the empirical value of the chiral condensate (qq) « —(250 MeV) 3 in the chiral limit [27], we choose n = 2 
for Eq. (11). This choice gives (qq) ~ —(239 MeV) 3 from Eq. (14) numerically which is comparable to the empirical 



value. Throughout the present work, we will use n = 2 for all the numerical calculations. 

Taking into account all the ingredients discusses so far, we arrive at the following neat expressions for the shear 
viscosity: 
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B. Temperature-dependent quark mass 



Here we would like to address how to determine the temperature dependence of the effective quark mass Mq in 
Eq. (11 ). In Refs. 31 , we derived it by using the caloron distribution with trivial holonomy, i.e. Harrington-Shepard 
29l [30] , Firstly, we want to explain briefly how to modify p and R as functions of T, using the caloron 
Details can be found in Ref. [3T]. An instanton distribution function for arbitrary N c and Nt can be written 



caloron 
solution 
with a Gaussian suppression factor as a function of T and an arbitrary instanton size p for pure-glue QCD 



d(p, T) = C Nc A b KS N ° p b ~ 5 exp [~(A N T 2 + ^np 2 )p 2 



(16) 



We note that the CP-invariant vacuum was taken into account in Eq. (16), and we assumed the same analytical 



form of the distribution function for both the instanton and anti-instanton. Note that the instanton number density 
(packing fraction) N/V = n = l/R 4 and p have been taken into account as functions of T implicitly. For simplicity, 
we take the numbers of the anti-instanton and instanton are the same, i.e. Nj = Nj = N. We also assigned the 
constant factor in the right-hand-side of the above equation as C for simplicity. The abbreviated notations are also 
given as: 
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Note that we defined the one-loop inverse charges ft and J3 at certain phenomenological cutoffs p cut and (R) ps -R. 
Ars denotes a scale, depending on a r enor malization scheme, whereas V3 for the three-dimensional volume. Using the 
instanton distribution function in Eq. ( 16 ), we can compute the average value of the instanton size p 2 straightforwardly 
as follows [40] : 
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where v = (b — 4)/2. It can be easily shown that Eq. (18) satisfies the following asymptotic behaviors [40] : 



lim p 2 (T) = ,/-= — , lim p 2 (T) = -. 



(18) 



(19) 



Here, the second relation of E q. (fl9| indicates a correct scale-temperature behavior at high T, i.e., 1/p ps A ex T. 
Substituting Eq. ( 18 ) into Eq. (116k, the caloron distribution function can be evaluated further: 



d{p, T)=C p b - 5 exp [- F{T)p 2 } , F{T) = \a n T 2 



AlT* 
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The instanton packing fraction n can be computed self-consistently, using the following equation: 

n^(T) = [Cr(i/)]', 



(20) 



(21) 



where we replaced NT/V3 — > n, and r(z/ ) stands for the T-function with an argument v. Note that C and /? can be 
determined easily using Eqs. (18) and ( |2l| ), incorporating the vacuum values for n w (200 MeV) 4 and p w (600 MeV) -1 : 
C ps 9.81 x 10~ 4 and /3 pa 9.19. Finally, in order for estimating the T-dependence of Mq, 
normalized distribution function, defined as follows, 



one needs to consider the 
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(22) 



Here, the subscript N denotes the normalized distribution. For brevity, we want to employ the large-A^ limit to 
simplify the expression for d^(p,T). In this limit, as understood from Eq. (22), diq{p,T) can be approximated as a 
(5-function: 



A', 



lim d N {p,T) = S[p-p{T)}. 



(23) 



The numerical result for the trajectory of p(T) is given in the panel (a) in Figurell] Here we choose p(0) = 1/A ps 1/3 
fm and R pa 1 fm for all the numerical calculations. These values are phenomenologically preferred in the present 
model [27]. The curve for p(T) shows that the average (anti)instanton size smoothly decreases with respect to 
temperature. It indicates that the instanton ensemble gets diluted and the nonperturbative effects via the quark- 
instanton interactions are diminished. At T = (150 ~ 200) MeV, which is close to the chiral phase transition 
temperature, the instanton size decreases by about (10 ~ 20)% in comparison to its value at T—0. Considering that 
the instanton size corresponds to the scale parameter of the model, i.e. UV cutoff mass, p ~ 1/A, the temperature- 
dependent cutoff mass is a clearly distinctive feature in comparison to other low-energy effective models, such as the 
NJL model. In addition, we also show the temperature dependence of the average (anti)instanton number density or 
(anti)instanton packing fraction, N/V, in the panel (a) of Figure]!] Again, the instanton number density decreases as 
temperature increases: The instanton ensemble is diluted. We will use these two temperature-dependent quantities 



for computing the shear viscosity in Eq. (15) 



C. Shear viscosity under strong external magnetic field 



Here, we briefly discuss how to induce the external magnetic field to the quark matter. Following the Schwinger 



method, we apply the minimal gauge substitution to the covariant derivative, id^ — > iD^ 



id u 



ie q Afj,. By doing 



that, the momentum dependent effective quark mass can be expanded in terms of the electric charge of the quark, 
then we obtain the following expression for 0(e q ) 39J: 
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We choose the specific configuration for the external magnetic field for convenience as 

B = (B x ,B y ,B z ) = (0, Bo sin 6 B , B cos B ), 



(24) 



(25) 



where 6b is an arbitrary angle. It has been verified that choosing arbitrary field configuration does not generate any 
qualitative difference. Considering that 1G = 1.95 x 10~ 14 MeV 2 in the natural unit and m 2 ps 10 18 G in terms of 
the pion mass m^ sw 140 MeV, it is quite convenient to employ the following parameterization for the magnetic field: 
eBo = nsm 2 . As for n — 1, the strength of the magnetic field is comparable to that of the magnetar. If ns becomes 
about (10 ~ 100), it can be compared to the strong magnetic field observed at the peripheral heavy-ion collisions at 
RHIC [25. 

Combining all ingredient we have the simple expression for the shear viscosity as a function of (T, Bq) up to 0(e 2 ): 
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where the summation runs over the light flavors u and d. Corresponding electrical quark charges are (e u ,e<i) — 
(+2/3, — l/3)e, in which e denotes the unit electrical charge e = V^rctEM i n the natural unit. Note t hat the magnetic 
field effect comes only from 3(e g i?o) 2 M| which is only proportional to Mq as understood as in Eq. (24). 



D. Calculating Entropy in our model 



As in Ref. [31 j . the LIM thermodynamic potential per volume in the leading 1/N C contributions at zero quark 
chemical potential can be written as follows: 
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where A and M represent a Lagrange multiplier to exponentiate the effective quark-instanton action and an arbitrary 
massive parameter to make the argument for the logarithm dimensionless. a stands for the isosinglet scalar meson 
field corresponding to the effective quark mass. In the leading 1/N C , we have the relation 2ct 2 = N/V [31]. The gap 
equation can be derived from Eq. (27) by differentiating f^uM by the Lagrange multiplier A: 
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Note that one can write the instanton packing fraction in terms of the effective quark mass M and p 



N 
V 



C N C M 2 

TT 2 /} 2 



(28) 



(29) 



The value of C locates in (1/3 - 1/4) for 1/p « 600 MeV, M « (300 ~ 400) MeV and N/V ra (200 - 260 MeV) 4 for 
vacuum gT]. We choose C = 0.27 to reproduce M = (340 — 350) MeV at (T, p) = in the chiral limit. After solving 
Eq. ( p8| with respect to Mo numerically, the numerical results for Mq as a function of T are given in the panel (b) of 
Figure |l| for the zero and finite current quark mass: m — (solid) and m = 5 MeV (dotted). These results indicate 
correct universal pattern for the phase transition pattern like the those of the Ising model, i.e. the second-order chiral 
phase transition for the massless fermion and the crossover for the finite mass. Here we choose the current quark 
mass to be about 5 MeV, considering the isospin symmetry for the light SU(2) flavor sector: m u sa m^ sa m = 5 MeV. 
From those numerical results, the phase transition temperatures for the two cases are obtained as To « (166, 170) 
MeV for 771 = (0, 5) MeV. The transition temperatures are indicated by the thin solid vertical lines in the panel (b) of 
Figurefl] Since we are interested in the ratio of the shear viscosity and the entropy density rj/s, we derive the entropy 
density s as follows: 
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FIG. 1: Average (anti)instanton size p ~ 1/A [fm] and (anti)instanton packing fraction (N/V) 1 ' 4 [GeV] as functions of T, 
computed from the Har ring ton- Shepard caloron distribution 29, 30 in the panel (a). Effective quark mass at zero virtuality, 
Mo computed from Eq. ( 28 1 as functions of T for m — (solid) and m = 5 MeV (dot), signaling the second-order and crossover 
chiral phase transitions, respectively, in the panel (b). The vertical lines indicate the chiral-phase-transition temperatures 
To = (166, 170) MeV for m = (0, 5) MeV. 



From the effective thermodynamic potential in Eq. (27 1, we obtain entropy density within the present model 
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In deriving Eq. (31 1, we assume that 2tr 2 sa N/V and AM ss A 4 as in the leading 1/N C) since A is only the scale 
parameter of the present model. The logarithm term ln[- • • ] in the first square bracket in the right-hand-side of Eq. ( 31 ) 



gives small contribution to the entropy density. As understood in the panel (a) of Figured] the (anti)instanton number 
density N/V is a function of temperature, so that its derivative with respect to T in the first them in right-hand-side 
of Eq. (31 1 is finite in general within the present model. The detailed calculations for these quantities will be given 



in a separated work 



III. NUMERICAL RESULTS AND DISCUSSIONS 



In this Section, we present and discuss the numerical results of the shear viscosity. In Figure [2j the shear viscosity 
is presented as a functions of T under the external magnetic field B = rigm^ in the chiral limit (a) and in the case 
of finite current-quark mass m — 5 MeV (b). The thick and thin lines indicate those with T-dependent parameters 
(TDP) and T- independent parameters (TIP), respectively. Note that TDP stands for that the scale parameter p and R 
are functions of T as shown in the left panel of Figure [l] whereas TIP means these two parameters are T-independent, 
i.e. p = p(0) and the same for R. The vertical lines shown in the both panels denote the transition temperatures T . 

We observe that the viscosity starts from zero, then keeps increasing as T increases for TDP, whereas it decreases 
beyond To for TIP. This observation suggests that the nontrivial T dependence of the scale parameter of a nonpertur- 
bative model makes significant effects on the behavior of the shear viscosity. As T becomes higher, the scale parameter 
p decreases implying the reduction of the SBxS effect in our theoretical framework: The (anti)instanton ensemble is 
diluted as T increases. In other words, the system is no longer strongly coupled one. It results in the larger shear 
viscosity. It is worth noting that the similar increasing behavior was also observed in the N JL model calculation [IT] . 
Although they considered small quark chemical potential n = 10 MeV and they treat the finite width for the quark 
spectral function as a free parameter. 

Furthermore, the shear viscosity becomes smaller under the strong magnetic field for the both cases of m = (0, 5) 
MeV. Again, this tendency can be understood by the enhancement of SB^S, in terms of the magnetic catalysis 143) . 
In our theoretical framework, the magnetic field contributions are proportional to M k ex Mq(T) as shown in Eqs. (24) 
and (126 1 . They disappear beyond T for the chiral limit as shown in the panel (a) of Figure pi due to the second-order 
chiral phase transition as understood by seeing the panel (b) in Figure nl the magnetic-field effects remain finite even 
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FIG. 2: (a) Shear viscosities r\ in the chiral limit as functions of T for different strengths for the static external magnetic 
field eBo — nsm^ for ub — (solid), 50 (dotted), and 100 (dashed) with the T-dependent parameters (TDP, thick) and 
T- independent parameters (TIP, thin), (b) The same curves with m = 5 MeV. The vertical lines indicate the chiral-phase- 
transition temperatures To — (166, 170) MeV for the (left, right) panels. 



beyond T in the case of the finite current-quark mass as in the panel (b) of Figure [2j It is because of the crossover 
pattern of the chiral restoration. At very high temperatures such as T > 220 MeV, the magnetic-field effects almost 
vanish even in the finite current-quark mass. Near the transition temperature Tq ~ 170 MeV, the shear viscosity 
becomes approximately r\ sa 0.02 GeV for all the cases. 

In literatures, the ratio of the shear viscosity and the entropy density rj/s has been considered as an important 
physical quantity. Hence we also present our result for rj/s here. First, in the left panel of Figure |3j we depict the 
entropy density using Eq. (31) for TDP (solid) and TIP (dash). Since we are interested only in the cases with the 
finite current-quark mass, we choose m — 5 MeV, as mentioned above. We find that rj/s are smoothly increasing 
curves with respect to T for the both cases, whereas the result for TDP is larger than the other. This can be easily 



understood by that the first term in the right-hand-side of Eq. (31 1 becomes zero for N/V = 1/i? 4 = const.. Note that 



here we set the external magnetic field to be zero here, since we have verified that the magnetic-field contributions to 
the entropy density are negligible. 

In the right panel of Figure pi we show the numerical results for the ratio rj/s as functions of T for TDP (thick) 
and TIP (thin), with different strengths of the magnetic field. Even the T dependence of the model parameters have 
been taken into account, the present model scale is about A f=a 600 MeV since it corresponds to the nonperturbative 
QCD region. Therefore, we confine our discussion to the temperature not much far beyond the chiral transition, i.e. 
T max = 350 MeV. The magnetic field dependence of rj/s comes only from the numerator rj. Again, here we choose 
the current-quark mass as m = 5 MeV which manifests the crossover chiral restoration pattern. The horizontal and 
vertical lines stand for the chiral transition temperature To = 170 MeV and the KSS-bound value rj/s — 1/(Att) « 0.08, 
respectively. The curves of rj/s for TDP decrease smoothly and approaches to the KSS bound. Those for TIP behave 
similarly but decrease more stiffly with respect to T. Note that the TIP curves undershoot the KSS bound at 
T f=a 270 MeV. These observations of the different curve behaviors implies that it is necessary to take the temperature 
dependence of the model parameters into consideration. The effects of the magnetic field are sizable below the chiral 
transition, then become negligible beyond T . Near the transition point, we observe only a few percent changes in 
the ratio rj/s due to the magnetic field. 

In the right panel of FigurepH the other theoretical estimations for the ratio rj/s is also presented for comparison. In 
Ref. [2TJ, the Monte-Carlo calculations of the two-point correlations in the pure SU(3) gauge were employed to compute 
the ratio with the nonperturbatively normalized operators. It gives rj/s — (0.134, 0.102) at T = (1.65, 1.24) T c , where 
T c is the critical temperature. This result is represented by solid square. The TDP curves are well compatible with 
their result at T sa 335 MeV, while the TIP curves undershoot the value. 

The effective models such as the NJL model have been also used for estimating the ratio as well. In Ref. [T7], it 
has been reported that rj/ ' s ~ (1/Att ~ 0.9) at (T, /j) = (200, 10) MeV, depending on the finite width for the quark 
spectral function. Averaging their values over the finite width, we have r\/ s sw 0.25, and it is represented by the solid 
circle in the left panel of Figure [3] It locates between the TDP and TIP curves. In their previous work with the 
same theoretical framework |llj . the shear viscosity increases with the larger quark chemical potential with a small 
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FIG. 3: (Color online) (a) Entropy density s as a function of T for with the T-dependent parameters (TDP, thick) and T- 
independent parameters (TIP, thin), (b) The ratio of the shear viscosity and entropy density rj/s in the same manner with 
the left panel, with different strengths of the external magnetic fields, ns — (0,50, 100), given in the (solid, dash, dot-dash) 
lines. We also showed the theoretical results from Meyer (LQCD) [2T] (square), Iwasaki (NJL) [T7] (circle), Sasaki (NJL) [19] 
(triangle), and Chen (xPT) [12] (diamond). The parameterization of the TDP curve for ub = in Eq. (32 1 is also given 
with the solid nabla. Detailed explanations for these theoretical values are given in the text. The vertical lines indicate the 
chiral-phase-transition temperatures To = 170 MeV for the (left, right) panels, while the horizontal one in the right panel stands 
for the lower bound of the QGP shear viscosity, i.e. the KSS bound x]j ' s — l/(47r) w 0.08. 



increasing rate. Hence the depicted point in the right panel is supposed to be lowered at ji—0. Nevertheless, the 
change from ji = (10 — > 0) MeV will not substantial in the present discussions. 

Employing the NJL model, Ref. [19] explored the transport coefficients near the chiral phase transition. From 
their results of the ratio rj/s « 0.5 at T sw 170 MeV. This value is depicted in the right panel of Figure [3] with the 
solid triangle. It is comparable with the TIP curves and but larger compared with the TDP curves by about 50%. 
Note that the rj/s curves in Ref. [19j show similar T-dependency to ours. However, beyond To, their curves turn into 
a slightly increasing ones and no longer similar to ours. In Ref. [12] they computed rj/s by using xPT below the 
chiral transition temperature. They estimated rj/s as a decreasing function of T with a typical value rj/s — 0.6 at 
T = 120 MeV with 50% uncertainty. We depict this value with the solid diamond with the error bar in the right 
panel of Figure [31 well-matching with the TDP curves. There are other theoretical estimations for rj/s for the high-T 
(T > 450 MeV) regions from LQCD and pQCD [331 HH HU [45] , and those results can not be reproduced within the 
present model: The theoretical results are usually larger than ours by (5 ~ 10) times. However, this situation is 
rather natural, considering that the present model is well compatible in the low-energy regions as a nonperturbative 
effective model. Hence, we do not discuss those result here. 

Finally, from the numerical results obtained above, we want to provide a simple parameterization of the ratio rj/s as 
a function of T. Since many theoretical approaches for the QGP dynamics have used a T- independent rj/s value [I], 
this parameterization would help to construct more realistic models of QGP. Taking into account that the magnetic- 
field effects are negligible for T > To as shown in the right panel of Figure [3j we just parameterize the numerical 
result for ns = 0. Employing a simple analytic form, one is led to 



T = (100 -350 MeV), 



(32) 



tj _ 0.87 1.19 0.28 

- =0.27- — + — -3-, 

s t t A t A 

where we use a notation t = T/T with T = 170 MeV. In the right panel of Figure pi we show the rj/s curve using 
the parameterization in Eq. (17) with the solid nabla. 



IV. SUMMARY, CONCLUSION, AND FUTURE PERSPECTIVES 



In summary, we have investigated the QGP shear viscosity of the SU(2) light-flavor quark matter at finite temper- 
ature under the strong external magnetic field. We employed the liquid instanton model and Green-Kubo formula to 
derive our result. The external magnetic field has been induced by the Schwinger method. Since the shear viscosity 
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becomes zero at the mean-field level, we suggested a finite-width quark spectral function motivated by the instanton 
model. There is one free parameter in this quark spectral function which has been determined by the chiral condensate 
value at zero temperature. The model parameters such as the average instanton size and inter-instanton distance are 
all temperature-dependent. This character is different from usual local-interaction models. We list our important 
observations as follows: 

• In our model, many parameters such as quark mass and instanton size are modified by temperature. Our way of 
modifying those parameters are supported by the fact that our effective thermodynamic potential at the leading 
large-iVc contributions manifests the correct chiral restoration patterns, i.e. the second order and crossover 
phase transitions for m = and m ^ 0, respectively. 

• Our study here has shown that the external magnetic field reduces rj due to the magnetic catalysis, i.e. the 
quarks are coupled more strongly in the presence of the magnetic field. The effect from the external magnetic 
field turns out to be sizable below the chiral transition temperature To = (166, 170) MeV for m = (0, 5) MeV. 
However it becomes inappreciable when temperature goes beyond T since the nonperturbative effects, such as 
the magnetic catalysis, becomes diminished. We also find that a typical value for the shear viscosity near To is 
r; = 0.02GeV 3 . 

• We observe that the T-dependent parameters (TDP), p(T) and R{T) play an important role beyond To which 
make 77 to increase. In contrast, rj starts to decrease after T , if we choose the T-independent parameters (TIP). 
The ratio of the shear viscosity and the entropy density, r\j s is also computed for the finite current-quark mass. 
It has been shown to be a monotonically decreasing function of T = (100 ~ 350) MeV. Furthermore, we also 
find that rj/s undershoots the KSS bound, rj/s = l/(47r), for TIP. On the other hand, rj/s approaches to the 
KSS bound for TDP. At T = 170 MeV, we find a typical value for the ratio as r\j s — 0.29 within the present 
model. 

• Our numerical results of r]/s for TDP are well comparable with other theoretical estimations such as the NJL 
model, LQCD, and x^T for T = (100 ~ 350) MeV, although we can not reproduce the very high-T results from 
LQCD and pQCD. However, this is rather natural, since the present model is well applicable for the low-energy 
regions. As for the future usage in QGP-dynamics studies, we also parameterize the numerical result of r]/s in 
a simple polynomial form as a function of t = T/Tq without the magnetic field. 

Encouraged by the fact our results obtained in the present work agree well with the empirical data, we would like 
to extend our study to other QGP transport coefficients, such as the bulk viscosity and heat conductivity. Moreover, 
it would be interesting to take into account the external electric field, which turns out to be considerably strong in 
heavy-ion collisions. Thus, the external electric field may bring considerable changes in the transport coefficients, in 
comparison to the pure magnetic one. Related works are under progress and will appear elsewhere. 
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Appendix 



The quark spectral function is normalized as follows: 

1 f x , . „ , f°° sgn(w) (W70 - a) 



2tt 



PFw(w,k)dw - 

sgn(w) (w7q - a) 
2\/2kEA 



exp 



2^2ttEA 
(w - Ef 
2A 2 



exp 



exp 



(w-E) 
2A2 

{w + Ef 



2 1 



2A 2 



exp 



(w + E) 
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2A 2 



dw 



dw. 



(33) 



Replacing the integral variable as w ± E = w± , Eq. ( 33 ) is led to 

,,2 



sgn(w + - E) [{w+ - E)"f Q - a] 
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exp 



2A 2 



sgn(w_ + E) [(w- + £)7o - a] 
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2K 2 



dw 
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sgn(-E)(-E"/Q -a) sgn(E){Ej - a) -{-Ej -a) +(£70 - a) 



2E 2E IE 2E 

which satisfies the spectral function normalization condition. 



70, (34) 
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